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Abstract 



A measurement of the Kg lifetime is presented using data recorded by the NA48 
experiment at the CERN-SPS during 1998 and 1999. The Kg lifetime is derived from 
the ratio of decay time distributions in simultaneous, collinear Kg and Kl beams, 
giving a result which is approximately independent of the detector acceptance and 
with reduced systematic errors. The result obtained is rs = (0.89598 ± 0.00048 ± 
0.00051) x 10" 10 s, where the first error is statistical and the second systematic. 
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1 Introduction 

Precise measurements of the basic physics parameters defining the neutral kaon 
system, such as the masses and mean lifetimes of the Ks and Kl states, are important not 
only in their own right but also as essential inputs to many kaon physics analyses such as 
studies of indirect and direct CP violation or precision tests of CPT invariance. 

The Ks lifetime is presently known with a relative precision of about 0.1% |]J, 
dominated by measurements from the NA31 experiment at CERN f|| and from the E731 
and E773 experiments at Fermilab [[], M. Here we present a measurement of the Kg lifetime 
from the NA48 experiment at the CERN-SPS, based on the same K° — ► data samples 
as used for the precise determination of the direct CP violation parameter Re(e'/e) ||. 
The Kg lifetime is measured using an analysis technique introduced by NA31, namely a 
fit to the Ks/Kl ratio of decay time distributions reconstructed in nearly collinear Kg and 
Kl beams. This gives a result which is essentially independent of the detector acceptance 
and therefore with reduced systematic errors; the Kl beam is in effect used to determine 
the detection efficiency for Kg decays. 

The analysis method is described in more detail in the next section, followed in 
Sections |3| and ^ by a brief summary of the main features of the NA48 experiment and of 
the reconstruction and selection of K° — > tctt decays; further details can be found in ||. 
The fit used to extract the Kg lifetime is described in Section |5] and the estimation of the 
systematic errors on the fitted lifetime is considered in Section |6|. 

2 The Method 

A defining principle of the NA48 experiment is the simultaneous recording of decays 
occurring within a common decay region traversed by two almost collinear beams of 
neutral kaons. The relative target positions for each beam, one far upstream of the decay 
region and the other much closer, ensure that kaon decays in the two beams are due 
dominantly to the Kl or to the Kg component, respectively. Assuming equal detection 
efficiencies for decays from the K s and K L beams, the ratio R = N S /N L of decay rates 
observed in each beam as a function of the longitudinal position z can be expressed as 

R(E K ,z) = A(E K )^\ (1) 

where is the kaon energy, A(E^) is a normalisation function which depends on the 
relative beam intensities and £s,l — ( z ~ -ZtargS^X^K/PK:) are the proper lifetimes for kaon 
decays in the Kg or Kl beams. The functions f(t) are given by 

/ SL (t) = e ~ t/Ts + M V /TL + 2D StL {E K )\r]\e- {t/TS+t/T ^ /2 cos (Am -t-<f>) (2) 

where rg and tl are the Kg and Kl mean lifetimes, Am = — mg is the mass dif- 
ference between the Kg and Kl states, \rj\ and <fi are the modulus and phase of the 
ratio A(Kl — > nn)/A(Ks — > nn) of decay amplitudes, and the dilutions Ds,l(Ek) = 
[JV(K°) - N(K°)]/[N(K°) + N(K° )] reflect the initial admixture of K° and K° in each 
beam. In practice, in place of z, it is convenient to analyse the longitudinal distribution of 
decay vertices in terms of the variable r = (z — zaksX^k/pk) which measures the decay 
proper lifetime relative to the upstream edge of the decay region, defined by the position, 
^aks; °f a set of scintillation counters. 

Except for the Kl beam at higher energies, where the interference term contributes 
appreciably, the functions /s(£s) an d are dominated by the exponential terms e~ is / TS 
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and 1 77 1 2 e W 71 , respectively. Hence the Ks/Kl ratio is approximately of the form 

R oc —— oc e-^V-s-i/n.) 
g— Wtl 

Since tl ^> rs, the ratio R is primarily sensitive to the Kg lifetime rg. The Kg lifetime 
is determined by fitting a function R(Ek,t) of the form given in Equations (|l|) and (§) 
to the ratio Ns/N^ of K° — > 7r7r decays reconstructed in the two beams. Besides the 
Kg lifetime itself, the normalisation A(.E'k) and the dilutions -Ds.iX-Ek) are also taken as 
free parameters in the fit, while the remaining physics parameters in equation (|2|) are 
taken from published measurements. Small acceptance differences between Kg and Kl 
decays are corrected using Monte Carlo, and background to the K L samples is subtracted 
using the data. The 7r + 7r~ and 7r°7r° decay modes are analysed separately and the results 
subsequently combined. 

The Kg/KL ratio is reconstructed in 20 bins of energy of width = 5GeV 

covering the range 70 < < 170 GeV and in lifetime bins of width At = O.lrg. The fit 
to determine the Kg lifetime is carried out in the lifetime range 0.5rg < r < 3.5rg. The 
lower lifetime limit of 0.5rg largely avoids detector resolution effects associated with the 
start of the decay region at r = while the upper limit at 3.5rg is dictated by the trigger 
requirements. The choice of lifetime range approximately minimises the total error on the 
measured Kg lifetime. 

3 The NA48 experiment 

The Kl and Kg beams in the NA48 experiment are derived from 450 GeV protons 
incident on separate targets positioned 126 m and 6 m upstream of the decay region, 
respectively || . The Kg target is located 72 mm above the axis of the Kl beam. The Kg 
target and collimator system is aligned along an axis directed slightly downwards such 
that the Kg and Kl beams subtend an angle of 0.6mrad and intersect at the centre of 
the detector, 120 m downstream of the Kg target. 

The beginning of the decay volume is accurately defined by a Kg anti-counter (AKS) 
located at the exit of the Kg collimator. The AKS is composed of a photon convertor 
followed by three scintillator counters and is used to veto all upstream decays in the Kg 
beam. To minimise interactions of beam particles with air and material, the decay region 
itself is contained in a 90 m long evacuated tank terminated by a thin composite polyamide 
(Kevlar) window of 3 x 10 -3 Xq thickness. The tank is followed by the NA48 detector, the 
principal components of which are a magnetic spectrometer for charged particle detection, 
a liquid krypton calorimeter for photon and electron detection, an iron-scintillator hadron 
calorimeter, and muon counters consisting of three planes of scintillator shielded by 80 cm 
thick iron walls. 

The charged particle spectrometer J7| consists of two drift chambers (DCH1 and 
DCH2) located before, and two drift chambers (DCH3 and DCH4) located after, a cen- 
tral dipole magnet. Each chamber has an area of 4.5 m 2 and is made up of four sets 
of two staggered sense wire planes oriented along four directions (horizontal, vertical, 
±45°). Track positions are reconstructed with a precision of 100 /im per view, while the 
momentum resolution is a(p)/p = 0.48% © 0.009 x p(GeV/c)%. 

The liquid krypton calorimeter (LKr) || consists of ~ 13000 cells, each of cross 
section about 2 cm x 2 cm and depth 27 radiation lengths. For the average photon energy 
of 25 GeV, the transverse spatial resolution for photon showers is 1 mm and the energy 
resolution is 0.8%. The energy response is linear to about 0.1% in the range 5—100 GeV. 
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The photon time resolution is ~500 ps and the 7r°7r° event time is known with a precision 
of ~220 ps. 

An evacuated beam pipe of radius 8 cm running along the full length of the detector 
on its central axis transports undecayed beam particles through each of the detector 
components. The beam tube and DCH drift chambers are aligned along the bisector of 
the converging Ks and Kl beams in order to equalise the acceptance for Kg and Kl 
decays. Undecayed neutral kaons in the Kl beam are largely confined within a transverse 
beam profile of radius ~ 3.5 cm at the position of the NA48 detector. The Kg beam has 
a larger beam divergence and a correspondingly larger transverse profile of about 5 cm 
radius. Due to scattering in the beam collimators, both the Ks and Kl beams have an 
associated halo of particles extending to larger radii from the beam axis. 

The trigger for 7r°7r° decays I5J requires that the total energy deposited in the LKr 
calorimeter be greater than 50 GeV, that the kaon impact point at the calorimeter (had 
it not decayed) be within 15 cm of the beam axis, that the decay vertex be less than 
5 Ks lifetimes from the beginning of the decay volume, and that each of the horizontal 
and vertical projections of the LKr energy distribution contain at most five peaks. The 

trigger operates with negligible deadtime and high efficiency, (99.920 ± 0.009)%, 
with no significant difference between Kg and Kl decays. The first level of the tt + tt~ 



trigger [ID] is based on signals from a scintillator hodoscope positioned in front of the 
LKr calorimeter and on the hit multiplicity in DCH1, and also requires a total energy in 
the LKr and hadron calorimeters of at least 35 GeV. The second level trigger is based on 
tracks reconstructed using the information from DCH1, 2 and 4, and includes requirements 
that the decay vertex be less than 4.5 Ks lifetimes from the beginning of the decay volume 
and that the reconstructed mass be larger than 0.95 7«k. The efficiency of the first and 
second level triggers is about 99.5% and 98.3%, respectively, again with no significant 
difference between K s and K L decays. The 7r + 7r~ trigger introduces a deadtime of about 
1.1%. 



4 Event Selection 

K° — > 7T7T decays are reconstructed and selected using the same procedures and 
selection requirements as for the Re(e'/e) analysis ||. The level of background remaining 
in the selected K° — > ttit samples, and the correction for acceptance differences between 
K s and K L decays, are also evaluated using similar techniques to those in ||. For the 
determination of the Kg lifetime, it is the lifetime dependence of these corrections which 
is of importance, rather than their energy dependence or their overall normalisation. 

4.1 The 7r°7r° sample 

The reconstruction of 7r°7r° events is based entirely on data from the LKr calorime- 
ter. Any group of four showers, each reconstructed within 5 ns of their average time, is 
considered. The energy of each shower is required to lie between 3 and 100 GeV, and show- 
ers close to the edges of the calorimeter (within 11 cm of the outer edge or within 15 cm of 
the central axis) or within 2 cm of a defective cell are excluded. The transverse distance 
between any pair of showers must be greater than 10 cm. The position of the centre of 
gravity of the event is defined as the energy-weighted average of the four shower positions. 
The radial distance, C g , of the centre of gravity from the detector axis is required to be 
less than 10 cm to suppress events due to the decay of particles in the beam halo. 

The kaon energy is estimated simply as the sum of the four shower energies. The 
longitudinal decay vertex position is reconstructed from the energies and positions of the 
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four showers, under the assumption that they come from the decay of a particle with the 
kaon mass, tok, moving along the beam axis. A resolution of ~ 40-60 cm is achieved on 
the vertex position, depending on the kaon energy and decay point. The shower pairing 
which best represents two 7r° — > 77 decays is inferred using a x 2 variable constructed 
from the sum, mi + m 2 , and difference, mi — m 2 , of the two candidate m 77 masses and a 
parameterisation of the resolutions on mi ± m 2 . 

Background from Kl — > 7r 7r°7r° decays is suppressed by requiring no additional 
showers with energy greater than 1.5 GeV within ±3ns around the event time, and by 
requiring \ 2 < 13.5 (which corresponds to 3.7a on the m 77 resolution). All other potential 
sources of background to the 7r°7r° sample are negligibly small. 

The level of K L — > 7r°7r 7r° background remaining in the K L — > 7r°7r° sample is 
estimated using a control region at large values of x 2 ■, 36 < x 2 < 135, dominated by 
background. Genuine Kl — > 7r°7r° signal events populating this control region are first 
subtracted using the Ks — > 7r 7r data, for which background is negligible. Small differences 
in the shape of the x 2 distribution for Ks — > 7r°7r° and Kl — > 7r°7r° decays are taken 
into account by applying a correction derived from Monte Carlo. The K L — > 7r°7r°7r 
background in the signal region, x 2 < 13.5, is estimated by extrapolating uniformly from 
the control region, with an additional factor A cxt = 1.2 ± 0.2 derived from Monte Carlo 
simulation. The background estimation is carried out separately for each bin of kaon 
energy and proper lifetime. At low energy and at low lifetime, the fraction of 7r 7r°7r 
background in the 7r°7r° sample is negligibly small, but rises with both energy and lifetime, 
reaching about 3% for E K = 170 GeV and r = 3.5r s . An increase of the background 
fraction with lifetime is to be expected since, when only four photons are detected in the 
LKr calorimeter, the reconstructed vertex position for a Kl — > 7r°7r 7r° decay is shifted 
downstream from its true position, while Kl — > 7r°7r 7r° decays occuring upstream of the 
decay region are largely removed by the K L beam collimators. 

4.2 The 7r + 7r _ sample 

K° — > 7r + 7r~ events are reconstructed from oppositely charged pairs of tracks found 
in the magnetic spectrometer. Each track is required to have momentum greater than 
10 GeV and to lie at least 12 cm from the centre of each DCH. Each track is also required to 
lie within the acceptance of the LKr calorimeter and the muon counters after extrapolation 
downstream. The momentum-weighted average of the track positions after extrapolation 
to the LKr calorimeter is used to define the position of the centre of gravity of the event. 
Its radial distance, C g , from the detector axis is required to be less than 10 cm. 

The separation between the two tracks at their point of closest approach after 
extrapolation upstream from the spectrometer is required to be less than 3 cm. The point 
midway between the tracks at their closest approach is used to define the decay vertex 
position. The resolution on the longitudinal position of the decay vertex is in the range ~ 
30-50 cm, while the transverse resolution is about 2 mm. The kaon momentum is computed 
from the opening angle 9 of the two tracks and from the ratio p+/p- of their momenta, 
assuming that the event corresponds to a K° — > it + it~ decay. Thus the reconstruction of 
both the decay vertex position and the kaon momentum, and hence also of the proper 
time for the decay, depend only on the geometry of the detector. 

Background to the 7r + 7r~ sample from A — > pn~ decays is reduced to a negligible 
level by applying an energy dependent upper cut on the track momentum asymmetry 
\p + —p-\/(p + + £>-)■ This cut also serves to remove events with tracks at low radius to 
the beam axis in a way which depends only on the momentum ratio of the two tracks. 
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Background from Kl — > 7rez/ (K e3 ) decays is suppressed by requiring E/p to be less than 
0.8 for each track, while Kl — > TTfiu (K^) decays are suppressed by rejecting events where 
one or both of the tracks is associated with a signal in the muon counters within ±4ns. 
Additional suppression of both K e3 and K^ 3 decays is obtained by requiring that the 7r + 7r~ 
invariant mass, m n7T , be compatible with the kaon mass to within 3cr of the (energy- 
dependent) mass resolution, and by requiring a small missing transverse momentum, 
pl^ < 200 MeV 2 /c 2 , where p' T is the component of the kaon momentum perpendicular to 
the line joining the production target (identified from the vertical position of the decay 
vertex, as described below) and the point where the kaon trajectory crosses the plane of 
the first drift chamber. The quantity has approximately the same resolution for 7r + 7r~ 
decays in both the K s and K L beams. 

The level of background remaining in the 7r + 7r~ sample is estimated by extrapolating 
an exponential fitted to the p^ distribution in a control region, 800 < p'^ < 2000 MeV 2 /c 2 , 
dominated by background, into the signal region p^ 2 < 200MeV 2 /c 2 . Due to limited 
statistics, the lifetime bin width was increased from O.Irs to 0.5rs for this procedure. 
The fraction of background in the Kl — > 7r + 7r~ sample is found to be about 2 x 10 -3 , 
and no significant dependence of the background fraction on lifetime is observed in any 
bin of energy. An appreciable lifetime dependence is not expected in this case since both 
the signal and background are observed in the detector as two-track final states with 
configurations which vary in similar fashion with the decay vertex longitudinal position. 



4.3 K s tagging 

For 7r + 7r~ decays, the good resolution on the transverse position of the decay vertex 
allows a clean separation of decays from the Kg and Kl beams; decays with a vertex 
position more than 4 cm above the Kl beam axis after extrapolation of the reconstructed 
parent kaon trajectory back to the position of the AKS counter were classed as belonging 
to the Kg beam ("vertex tagging"). For 7r°7r° decays, only the longitudinal position of 
the decay vertex can be reconstructed. In this case, the identification of decays from 
the Ks beam is accomplished using a tagging station (Tagger) traversed by the proton 
beam during transport to the Ks target and consisting of two scintillator ladders (one 
horizontal, one vertical) fi~l |. Each scintillator counter has a time resolution of ~140ps, 
and a proton crosses at least two counters. A 7r°7r° decay is classed as belonging to the 
Ks beam if a signal is observed in the Tagger within ±2ns of the reconstructed event 
time. A small fraction, = (1.6 ± 0.5) x 10~ 4 , of Ks — > 7r°7r decays are mistagged 
as belonging to the Kl beam due to Tagger detection inefficiencies. A larger fraction, 
a LS = (10.692 ± 0.020)%, of Kl — > tt tt° decays are mistagged as belonging to the Ks 
beam due to accidental proton signals in the Tagger in time with the event. The values of 
and c^ls were inferred from studies of the mistagging fractions and a^s" for vertex 
tagged 7r + 7r~ decays, combined with studies of Kl — > and Kl — > decays to 

determine the small differences in mistagging rates between charged and neutral modes. 



4.4 Event samples 

Applying the selection cuts above to the data recorded by NA48 during 1998 and 
1999 yielded samples of 13.2M K s -> n+it~, 12.2M K L -> tt+tt-, 3.1M K s -> 7r°7r° and 
2.8M Kl — > 7r°7r° events with reconstructed lifetimes in the range 0.5 < r/Vs < 3.5, where 
the number of 7r°7r° events has been corrected for Ks-Kl mistagging. 

Samples of simulated K° — > tttt events corresponding to about 3-5 times the data 
statistics were also available for analysis. The Monte Carlo included a detailed modelling 
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of the Ks and Kl beams (including the Ks beam halo) and used the GEANT package |T2| 
for particle tracking and the simulation of processes such as multiple scattering, photon 
conversion, and secondary interactions. Simulation of the LKr calorimeter response was 
based on a library of electromagnetic and hadronic showers generated using GEANT. The 
simulated events were passed through the same reconstruction and selection code as the 
data. Asymmetric non-Gaussian tails in photon shower energies due to hadron photopro- 
duction in the liquid krypton, which arise in about 3 x 10~ 3 of cases, were modelled using 
a parameterisation which was applied to the reconstructed shower energies in the Monte 
Carlo events. 

4.5 Acceptance correction 

Since the Ks and Kl beams are almost collinear, and since Ks and Kl decays are 
recorded simultaneously using a common trigger, the acceptances for Ks and Kl decays 
as a function of energy and position are equal to good approximation, and the detector 
acceptance essentially cancels in the Ks/Kl ratio of lifetime distributions from the two 
beams. For the 7r + 7r~ decay mode, effects due to decay in flight of the charged pions also 
cancel in the Ks/Kl ratio. 

Small acceptance differences arise due to the different divergences and transverse 
profiles of the Ks and Kl beams and the consequent differences in illumination of the 
detector. Also, the different definition of the upstream edge of the decay region for the 
two beams introduces large acceptance differences at low values of the reconstructed 
lifetime. The AKS veto requirement effectively represents a cut r > on the true lifetime 
of decays in the Ks beam while the upstream acceptance for Kl decays extends smoothly 
into the region r < 0. This difference is clearly seen in Figure [I] where examples of the 
uncorrected reconstructed lifetime distributions for the Ks and Kl beams are plotted. 

Acceptance differences for Ks and Kl decays were accounted for by dividing the 
Ks/Kl ratio observed in the data by the acceptance ratio predicted by the Monte Carlo. 
The acceptance correction also took detector resolution effects (smearing) into account 
by defining the Ks and Kl acceptances as the ratio of the number of simulated events 
passing the selection cuts in a bin of reconstructed energy and lifetime to the number of 
generated events in the corresponding bin defined using the true energy and lifetime. 

Resolution induced acceptance differences due to the AKS veto requirement are 
largely eliminated by restricting the lifetime range considered to the region r > 0.5rs. A 
small (up to ~ 2%) residual variation of the acceptance correction remains for 7r°7r° events 
at low energy and lifetime, but the correction for 7r°7r° events is otherwise independent of 
lifetime within the errors due to Monte Carlo statistics. For -k + tt~ events, the Ks/Kl ac- 
ceptance ratio is independent of lifetime at lower energies, but falls significantly (by up to 
~ 7%) with increasing lifetime at higher energies. This is not due to an intrinsic difference 
between the Ks and Kl acceptances, but rather to a slow decrease of the acceptance with 
increasing radial distance of the decay vertex from the beam axis, convoluted with the 
larger transverse size of the Ks beam. The fall in acceptance with increasing radius is due 
dominantly to the geometrical cuts and to the cut on the track momentum asymmetry, 
and can be reliably modelled and checked by comparing the relevant reconstructed distri- 
bution in data and Monte Carlo. Similarly, the modelling of the transverse beam profiles 
is checked by comparing the reconstructed C g distributions in data and Monte Carlo. 
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5 Results 

The Ks/Kl ratio R = Ns/N^ is computed in two-dimensional bins of energy and 
lifetime of width AEk = 5GeV and At = O.Irs, and corrected bin-by-bin for residual 
background in the Kl — > 7r + 7r~ and Kl — > 7r°7r° samples and for small differences in the 
K s and K L acceptances, as described above. For the 7r°7r° mode, the Kg/K L ratio is also 
corrected for Ks-Kl mistagging. Examples of the dependence of the corrected Ks/Kl ratio 
on the reconstructed proper lifetime r are shown in Figure [I]. The shapes of the separate 
Ks and Kl lifetime distributions change appreciably with energy, due largely to changes 
in the lifetime dependences of the 7r + 7r~ and 7r°7r° acceptances, but the corrected Ks/Kl 
ratio is approximately independent of energy. 

The Ks lifetime is determined by fitting a function of the form given in Equations ([!]) 
and (0) to the Ks/Kl ratio over the lifetime range 0.5 < r/rs < 3.5 and energy range 
70 < Ek < 170 GeV. The following x 2 quantity is minimised in the fit: 



where £7j is the central energy of the z'th bin of energy, tj is the central value of the j'th bin 
of proper lifetime, R^ is the corrected Ks/Kl ratio from the data, is the statistical error 
on R^, and R(Ei,Tj) = A/s(£s)//l(£l) is the expected Ks/Kl ratio. The normalisation 
parameters A\ = A(Ei) are allowed to vary in the fit. At low energies, 70 < Ek < 140 GeV, 
the values of the dilutions D$(Ek) and Dl(Ek) appearing in the interference terms in 
the functions /s(£s) an d are fixed using measurements of the dilution by the NA31 

experiment |]13|| . The dilutions D^(Ek) and D^(Ek) rise from about 0.2 at Ek = 70 GeV 
to about 0.35 at Ek = 140 GeV and differ slightly because of the different kaon production 
angles (4.2mrad and 2.4mrad) used to define the Ks and Kl beams; small corrections are 
applied to the NA31 measurements to account for the different production angles used 
in this experiment. At high energies, 140 < Ek < 170 GeV, where the interference term 
contributes significantly to the ratio R(E k ,t), the values Dj = D s (Ei) = D L (Ei) of the 
dilution are allowed to vary in the fit, neglecting the small differences between the Ks 
and Kl dilutions. The remaining physics parameters appearing in the functions f(t) are 
taken from the PDG averages of existing published measurements, except the Ks lifetime 
Ts which is a free parameter in the fit. 

The fit was carried out separately for the 7r + 7r~ and 7r 7r samples from each year of 
data taking. The values of the Ks lifetime from each fit are summarised in Table [TJ. The 
errors are statistical only and include the contribution from finite Monte Carlo statistics. 
The minimum value of the fit \ 2 is a l so given in Table [I] and corresponds in all cases to an 
acceptable fit probability. Examples of the fit results and the fit quality for the and 
7r°7r° data samples from the 1999 run are shown in Figure |l|. The effect of the interference 
terms in the fit function is visible in the plots for the high energy bin as a slight curvature 
of the lines representing the fit result. 

Consistent values of the fitted Ks lifetime rs, and of the fitted normalisation and 
dilution parameters Ai and Di, were found for the charged and neutral modes and for 
the samples from the different years of data taking. The fitted values of the dilution Di 
within the energy range 140 < Ek < 170 GeV were also consistent with the published 
NA31 measurements. 
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Figure 1: The points with error bars in the upper plot of each pair show examples of the 
corrected Ks/Kl ratio as a function of the reconstructed proper lifetime, r, expressed in 
nominal K s lifetime units of 0.8927 x 10~ 10 s. The dotted and dashed histograms show 
the corresponding uncorrected lifetime distributions for the Ks and Kl beams, arbitrarily 
normalised. The curves show the results of the fit for the Ks lifetime, while the lower plot 
of each pair shows the normalised fit residuals. 
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7T + 7T 


7T 7r 


Tg/10- 1U S 


rVdof 


r s /10- lu s 


X 2 /dof 


1998 
1999 


0.89578 ±0.00109 
0.89598 ±0.00072 


628.2/573 
601.2/573 


0.89606 ± 0.00247 
0.89635 ±0.00167 


551.0/573 
543.5/573 



Table 1: Values of the Kg mean lifetime and the minimum chi-squared per degree 
of freedom from the lifetime fits. The errors are statistical only. 



6 Systematic errors 

Various sources of systematic error on the measured Kg lifetime were considered, 
including uncertainties in the reconstruction of the energy and distance scales for the 
charged and neutral decay modes, and uncertainties due to the background subtraction, 
the tagging correction, the Monte Carlo acceptance correction, and the external physics 
parameters used in the fitting function. 

As noted in Section fO|, the energy and distance scales for reconstructed 7r + 7r~ events 
are determined largely by the detector geometry, especially of DCH1 and DCH2. The 
systematic error due to uncertainties in detector geometry was estimated by considering 
a variation of ±2 mm in the longitudinal separation of DCH1 and DCH2 and a variation 
of ±20 / um/m in the relative transverse scale of the two chambers. The position of the 
upstream edge of the K s — > 7r + 7r~ decay vertex distribution was found to be consistent 
within these tolerances with the nominal position of the AKS counter. In addition, the 
momentum scale for reconstructed tracks was varied by ±0.1%, as determined from the 
consistency of the reconstructed 7r + 7r~ invariant mass with the kaon mass, wk. In each 
case, the kaon decay reconstruction was redone for all events in the data using the modified 
track positions or momenta, and the lifetime analysis repeated. The resulting change in 
the fitted Kg lifetime was taken as the estimate of the corresponding systematic error. 

The systematic error due to uncertainties in the reconstruction of 7r°7r° events was 
estimated by varying the energy scale of reconstructed showers by ±0.03% and by varying 
the linearity and uniformity of the calorimeter response by modifying the reconstructed 
shower energy by an amount AE = a + (3E 2 ± ^rE, where r is the radial distance of the 
shower from the central detector axis. The allowed variation in the neutral energy scale is 
determined from a comparison of the position of the edge of the Kg — > 7r°7r° decay vertex 
distribution with the nominal position of the AKS counter. The allowed ranges of the 
constants a, (3 and 7 are determined from studies of K e 3, K° — > 7r°7r° and Kl — > 7r 7r 7r° 
decays, and of 7r° — > 77 and r\ — * 77 decays in special runs with n~ beams. 

The uncertainty due to the Kl — > 7r + 7r~ background subtraction was assessed by 
conservatively assigning a systematic error equal to ±100% of the change in the fitted Kg 
lifetime when the background subtraction was removed. The systematic error associated 
with the Kl — > 7r°7r° background subtraction was estimated by varying the extrapolation 
factor A ex t within the range A ex t = 1.2 ± 0.2. 

The Kl — > 7r + 7r~ and Kl — > 7r°7r° samples contain a small fraction of events aris- 
ing from the genuine irir decay of neutral kaons produced by scattering in (dominantly) 
the final collimator of the Kl beam. In the charged mode, collimator scattered events 
are largely removed by the cut on p^ 2 and contribute a fraction less than 10~ 4 to the 
Kl — > 7r + 7r~ sample. This contribution is already included in the charged background es- 
timate described in Section |47|, and no additional systematic error is warranted. Studies 
of events with large transverse momentum (p^ > 200MeV 2 /c 2 ) showed that the lifetime 
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Source 


Variation 


7T + 7T 


7T 7r° 


7T7T 


DCH1 radial scale 


±20 /im/m 


±1.8 




±1.5 


DCH1 z position 


±2 mm 


±0.2 




±0.1 


Momentum scale 


±0.001 


±0.3 




±0.3 


LKr Energy scale 


±0.0003 




±1.7 


±0.3 


Energy loss (a) 


±10MeV 




±3.4 


±0.6 


Non-linearity (/3) 


±0.00002 GeV- 1 




±1.6 


±0.3 


Non-uniformity (7) 


±0.00001 cm" 1 




±1.0 


±0.2 


LKr radial scale 


±0.0003 




±3.0 


±0.5 


Charged background 


±100% 


±1.4 




±1.2 


Neutral background 


A ext = 1.2 ± 0.2 




±2.1 


±0.3 


Collimator scattering 


±100% 




±3.2 


±0.5 


MC: Kg beam y posn 


±2 mm 


±1.6 


±0.2 


±1.3 


MC: Kl beam y posn 


±2 mm 


±1.4 


±0.2 


±1.2 


MC: K s beam halo 


±100% 


±1.3 


±1.3 


±1.4 


MC: non-Gaussian tails 


±50% 




±4.4 


±0.7 


MC: statistics 




±3.0 


±4.9 


±2.7 


r L 


5.17 ±0.04 x 10" 8 s 


±0.1 


±0.1 


±0.1 


\v+-\, \Voo\ 


2.276, 2.262 ±0.017 x 10" 3 


±0.1 


±0.1 


±0.1 


0+-, 000 


43.3 ±0.5°, 43.2 ± 1.0° 


±0.8 


±0.4 


±0.8 


Am 


0.5300 ±0.0012 x lO^s- 1 


±0.5 


±0.7 


±0.5 


D(E K ) 


±2 a 


±0.3 


±1.0 


±0.4 




0.10692 ±0.00020 




±2.1 


±0.4 


"SL 


(1.6 ±0.5) x 10" 4 




±0.4 


±0.1 


Fit method 




±2.7 


±2.7 


±2.7 


Total sys 


:ematic error 


±5.4 


±10.9 


±5.1 


Statistical error 


±5.2 


±12.9 


±4.8 



Table 2: Summary of systematic errors on the measured Kg lifetime, in units of 10 14 s. 
The final column corresponds to a combination of the 7r + 7r~ and 7r°7r° results. 



distribution of collimator scattered events corresponds approximately to an exponential 
decay characterised by the Kg, rather than Kl, lifetime, and that their energy distribution 
is similar to that of unscattered K L decays. The high events were also used to estimate 



the fraction, (4.2±1.0) x 10 [141, of collimator scattered events in the Kl — > tt tt sample 



where no transverse momentum cut is possible. For the 7r°7r° mode, the effect of collimator 
scattering on the Kg lifetime measurement was studied by modifying the fit function to be 
of the form R(Ei, Tj) = Aif s (t s )/[f L (t L ) + A co ne" Wl/rs) ] where t co n = (z - z coU )(m K /p K ) 
is the proper lifetime of the decay relative to the position of the final Kl collimator and 
Aeon is a constant which is adjusted to give the above fraction of collimator scattered 
events in the Kl — > 7r°7r° sample. This fraction was assumed to be independent of energy. 
The resulting systematic error was conservatively estimated as ±100% of the change in 
the Kg lifetime due to inclusion of the collimator scattering component in the fit. 

The systematic error associated with uncertainties in the acceptance correction was 
estimated by varying the assumed vertical position of the Kg and Kl beams in the Monte 
Carlo within a range ±2 mm around their nominal positions. In addition, the sensitivity 
to the transverse beam profile was studied by switching off the simulation of the Kg beam 
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halo in the Monte Carlo. In each case, the resulting change in the fitted Ks lifetime was 
taken as the systematic error. The sensitivity of the acceptance correction to the modelling 
of the detector resolution was assessed by removing the simulation of non-Gaussian tails 
from the reconstructed photon shower energies, and a systematic error of ±50% of the 
effect on the fitted Ks lifetime was assigned. The effect of possible non-Gaussian tails in 
the reconstructed drift chamber hit positions was also studied in the Monte Carlo, but 
found to have a negligible effect on the fitted Ks lifetime. Finally, the component of the 
overall statistical error on the fitted lifetime arising from the finite Monte Carlo statistics 
was extracted, and classified as a separate systematic error. 

The external physics parameters tl, t)+-, t]oo, 0+-, 4>oo and Am appearing in the 
functions f(t) were each varied in turn within a range given by the error on the PDG 
average value. The mistagging fractions a^ s and «sl were varied within the uncertainties 
given in Section fO| . The value of the dilutions Ds,l(Ek) in the energy range 70 < Ek < 



140 GeV was varied by twice the error on the NA31 measurements J0|]. The extra factor of 
two conservatively takes into account uncertainties in correcting the NA31 measurements 
to the NA48 experiment. 

The fitting procedure itself was tested by applying the fit to the generated inclusive 
energy and lifetime distributions of parent kaons from the Monte Carlo samples. No 
significant bias was observed on the fitted Ks lifetime, and the statistical precision of the 
test was assigned as a systematic error. 

The systematic errors on the Ks lifetime from each of the above sources for the 
1998 and 1999 data combined are summarised in Table |2]. For comparison, the effect of 
the background subtraction is to increase the fitted Ks lifetime for the 7r°7r° mode by 
12.8 x 10 _14 s, while the acceptance correction from the Monte Carlo changes the Ks 
lifetime by +22.5 x 10~ 14 s and —9.1 x 10 -14 s for the ti + tt~ and 7r°7r° modes, respectively. 

The analyses carried out for each year of data taking and for each decay mode 
are combined taking into account any correlations between the separate analyses. The 
total systematic error is obtained by summing the individual errors in quadrature. The 
measured values of the Ks lifetime for each decay mode and for the two decay modes 
combined are: 

r s = (0.89592 ± 0.00052 ± 0.00054) x 10~ 10 s (tt+tt") 

10 „ /„0 0\ 



r s = (0.89626 ± 0.00129 ± 0.00100) x 10~ iU s ( 



7T 7T 



r s = (0.89598 ± 0.00048 ± 0.00051) x 10~ iU s (tt+tt" + ttV ) 

where the first error is statistical and the second systematic. 

Various cross-checks were performed to verify the integrity and stability of the re- 
sult. No significant dependence of the fitted Ks lifetime was found on the lifetime or 
energy range used in the fit, or on the energy range within which D(Ek) was allowed 
to vary in the fit. The effect of varying the main selection cuts used in the analysis was 
also studied, namely the cuts on p^, C g , the momentum asymmetry \p + —p~\/ (p+ +P-), 
and the minimum radii of tracks and clusters in DCH1 and the LKr calorimeter. In each 
case, either no statistically significant variation in the fitted lifetime was found, or the 
observed variation was found to be within the systematic errors assigned. The stability of 
the analysis was also tested by dividing the 7i + n~ data sample into topologies for which 
the positive and negative tracks initially curve towards or away from each other in the 
spectrometer magnetic field; consistent values of the fitted lifetime were found for the two 
topologies. Similarly, selecting either of the two polarities of the spectrometer magnetic 
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field setting gave consistent results. Other tests involved dividing the data samples ac- 
cording to the primary data-taking periods, the event time within the SPS spill, and the 
azimuthal orientation of the decay. For the charged mode, separating Kg and Kl decays 
using the Tagger in place of vertex tagging gave no significant change in the result. 

7 Summary 

The Kg lifetime has been measured using Kg — > and Kl — > decays recorded 
by the NA48 experiment in 1998 and 1999. The combined result for the 7r + 7r~ and 7r°7r° 
decay modes, r s = (0.89598 ± 0.00048 ± 0.00051) x 10" 10 s, has a precision better than 
that of the current PDG average JI| of existing measurements, and lies about 1.7 standard 
deviations above it. 
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